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Using the Kantorovitch method in combination with a Gaussian ansatz, we derive the equations of motion 
for spatial, temporal and spatiotemporal optical propagation in a dispersive Kerr medium with a general 
transverse and spectral gain profile. By rewriting the variational equations as differential equations for the 
temporal and spatial Gaussian q parameters, optical ABCD matrices for the Kerr effect, a general transverse 
gain profile and nonparabolic spectral gain filtering arc obtained. Further effects can easily be taken into 
account by adding the corresponding ABCD matrices. Applications include the temporal pulse dynamics in 
gain fibers and the beam propagation or spatiotemporal pulse evolution in bulk gain media. As an example, the 
steady-state spatiotemporal Gaussian pulse dynamics in a Kerr-lens mode-locked laser resonator is studied. 
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I. INTRODUCTION 



The optical propagation in Kerr media with a trans- 
verse and spectral gain profile is of much interest in many 
areas. For instance, a combination of gain guiding and 
nonlinear self-focusing can play an important role for 
the spatial propagation of a laser beam in high-power 
laser rods. In nonlinear dispersive gain fibers, e.g., in 
fiber lasers or optical transmission lines equipped with 
erbium-doped fiber amplifiers, the temporal pulse evo- 
lution is affected by self-phase modulation and spectral 
filtering. 1 The coupled spatiotemporal dynamics in the 
gain medium is relevant for the operation of Kerr-lens 
mode-locked (KLM) lasers, where the pulse stabilization 
is governed by a combination of spatial and temporal 
effects^— 

The variational approach has been extensively used for 
an approximate description of the optical propagation 
in nonlinear media. By describing the optical field in 
terms of a trial function with free parameters, a set of 
coupled ordinary differential equations can be extracted 
from the partial differential equation governing the op- 
tical propagation. This allows for an analytical analysis 
or an efficient numerical treatment using a standard dif- 
ferential equation solver. The Rayleigh-Ritz method is 
a widely-used variational technique for the treatment of 
conservative systems, and has been applied to the spa- 
tial, temporal and spatiotemporal optical propagation 
in Kerr media. 5-9 Different approaches have been de- 
veloped to include dissipative effectsJ^— Here, we use 
a generalization of the Rayleigh-Ritz method known as 
Kantorovitch method* 1 ^ This technique has for example 
been applied to the nonlinear temporal pulse propaga- 
tion including parabolic spectral gain filtering,— ^ using a 
Pereira-Stenflo type ansatz, and to Gaussian beam prop- 
agation in air. 14 

In the following, we apply the Kantorovitch method to 
the description of the full spatiotemporal optical prop- 
agation of Gaussian light bullets in Kerr media, taking 
into account an arbitrary gain profile. As temporal ef- 
fects, we consider dispersion, self-phase modulation and 



spectral gain filtering. The spatial effects include diffrac- 
tion, self-focusing and a transverse gain profile. Also the 
cases of purely spatial beam propagation and purely tem- 
poral pulse evolution are considered. The assumption of 
a Gaussian ansatz allows us to relate the equations of mo- 
tion to the compact and elegant ABCD matrix formal- 
ism for Gaussian beam and pulse propagation! 15 i 16 By 
rewriting the variational equations as differential equa- 
tions for the q parameters, we can extract ABCD matri- 
ces for the Kerr effect, a general transverse gain profile 
and nonparabolic spectral gain filtering. Further effects, 
like a parabolic refractive index profile, can easily be in- 
corporated by adding the corresponding matrix elements. 

The paper is organized as follows: In Section [Til the 
spatial, temporal and spatiotemporal equations of mo- 
tion for the Gaussian parameters are obtained from the 
generalized nonlinear Schrodinger equation, which gov- 
erns the optical propagation in Kerr media. In Section 
III! these equations are reformulated within the frame- 
work of the ABCD matrix formalism, taking advantage 
of the Gaussian q parameter description. As an example, 
the Gaussian pulse dynamics in a KLM laser resonator is 
studied in Section HVl including a soft gain aperture and 
spectral filtering. The paper is concluded in Section |Vl 



II. VARIATIONAL APPROACH 

A linearly polarized light pulse, propagating in z direc- 
tion through a dispersive Kerr medium with a parabolic 
transverse and spectral gain profile, is described by the 
generalized nonlinear Schrodinger equation 

id z U - Vd 2 U + B {dl + flg) U + 5 \U\ 2 U = Q (1) 

with the gain term 



<3 = i (so - 9xx 2 - g y y 2 + g u d^) U. 



(2) 



The retarded time is defined as t T = t — z/v s with the 
group velocity v g , x and y are the transverse coordinates. 
U is the slowly varying envelope, normalized such that 
its absolute square gives the intensity of the wave. The 
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transve rse electr ical field component E t is related to U by 
E t = ^2Zq/uq x 5R{t/exp [i (kz — wot)]}, where no and 
k = ngko are the refractive index and the wavenumber 
at the center frequency ujq, and Zq is the wave resis- 
tance in vacuum. Here we use the 'physics' convention, 
in which a plane wave is described by exp [i (kz — u)at)]. 
The 'engineering' notation exp [j (wot — kz)] can easily be 
obtained by the formal transcription i — > — j in all expres- 
sions. The parameters for dispersion and diffraction are 
given by T> = \k" , where k" is the second derivative of 
the wavenumber at loq, an d 6=1/ (2k). The nonlinear- 
ity parameter is 5 = k^nn, where n^i is the intensity 
dependent refractive index, so that the total refractive 
index is given by n = no + nil \ U\ . In general, the co- 
efficients V, B, S, go, g x , g v , gu, and thus no and k" ', 
depend on the position z in the medium. For example, 
the material parameters change abruptly at the interface 
between two materials, and in optically pumped gain me- 
dia, go, g x and g y depend on z due to the divergence of 
the pump beam. Although the z argument is suppressed 
for a more compact notation, all the equations given in 
this paper are valid for z dependent coefficients. 

As a trial function for U, we choose a complex Gaus- 
sian, which has been widely used for the variational anal- 
ysis of optical propagation in Kerr media^^— Since it is 
an exact solution of Eq. |T]) for S = 0, we expect it to 
be a good approximation to the exact solution, at least 
for moderate nonlinearity. It has been shown that for 
the temporal dispersion-managed soliton dynamics, the 
Gaussian description is applicable over a wide param- 
eter rangej 17 ' 18 The transverse Gaussian field distribu- 
tion, which is the fundamental mode in linear paraxial 
resonators, is widely used as an approximate descrip- 
tion for the transverse field distribution in nonlinear 
resonators i 3 ' 4 ' 19 In addition, the Gaussian trial function 
can conveniently be characterized in terms of the com- 
plex q parameters, which allow for a compact description 
of the optical propagation based on the ABCD matrix 
formalism i 15 ! 16 

In the following, the Gaussian equations of motion, 
obtained by the Kantorovitch method, are given for the 
spatial beam propagation and the temporal pulse evolu- 
tion in a Kerr medium with a parabolic transverse and 
spectral gain profile, as well as for the full spatiotem- 
poral dynamics. The derivation of the spatiotemporal 
equations can be found in Appendix A; the purely spa- 
tial and temporal equations are obtained in an analogous 
manner. 



A. Equations of Motion 

The generalized nonlinear Schrodinger equation, as 
given in Eq. fTJ, describes the spatiotemporal pulse dy- 
namics in a dispersive Kerr medium, for example the 
propagation of a light bullet in a Kerr-lens mode-locked 
(KLM) laser, taking into account a parabolic transverse 
and spectral gain profile. The Gaussian test function is 



given by 

U(z,t T ,x,y) 



U (z) exp 
1 

2«£(z) 
1 

2^Tz) 



2T 2 (z 
ia x (z) 



- [a y( z ) 



y 



(3) 



with the beam widths w x and w v , the pulse duration T, 
the spatial chirp parameters a x and a y , and the tempo- 
ral chirp parameter b. The complex amplitude can be 
written as 



U(z) = A(z) exp [i(/)(z)} 



(4) 



The derivation of the equations of motion is given in Ap- 
pendix A. The equations for the beam width and the 
pulse duration are given by 



w„ 



4Ba p Wp 



T' = -<WbT 



9u 



— - Ab 2 T 3 



(5a) 
(5b) 



where p = x,y, and the prime denotes a partial deriva- 
tive with respect to z. Taking the full spatiotemporal 
dynamics into account, the intensity dependent contri- 
butions in the equations for the chirp parameters and 
the phase have different prcf actors c a>r j,, as compared to 
the purely temporal or spatial dynamics. For the chirp 
parameters, we obtain 
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b' = —T> [ -^-r — 4b 2 

rp4 



A 2 

c a 5^r, (6a) 
w 

p 

b A 2 
-43 w — -c a< 5— , (6b) 



with p = x,y and c a = V2/8. For the amplitude, the 
variational principle yields 



A' 



A (go 



9lo 

T 2 



2Vb - 2Ba x - 2Ba, 



(7) 



and the phase evolution is described by 
1 



6' = 2gJ> + V 



rp 2 



bX-bX 



+ C4> SA 2 (8) 



with c = 7V2/16. 

Setting V = g u = in Eqs. (JXJ) and ^ yields the 
nonlinear Schodinger equation for the purely spatial dy- 
namics. This equation describes the cw propagation of a 
beam with a transverse field distribution U = U (z, x, y) 
in a Kerr medium with a transverse gain profile, for ex- 
ample the nonlinear gain medium of an optically pumped 
solid-state laser. The Gaussian trial function is given by 



U(z,x,y) 



U (z) exp 
1 

2w 2 (z) 



1 



2wl(z) 



ia y (z) 



(9) 
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with the complex amplitude defined in Eq. (QJ. The rele- 
vant equations of motion are here Eqs. (|5aj) , (|6al) . ([7]) and 
([8|) with T> = g u =0. The nonlinearity coefficients are 
now given by c a = 1/4 and = 3/4, which are the same 
as derived using the method of minimum weighted square 
mean error,— but different from the results obtained by 
a Taylor expansion^ 

The propagation equations for the purely temporal 
pulse dynamics are obtained by setting B = g x = g y = 
in Eqs. ^ and ([2]). This equation describes the propa- 
gation of a pulse U — U (z,t T ) in a nonlinear dispersive 
medium with frequency dependent loss or gain, like an 
optical amplifier. This equation is also referred to as 
complex cubic Ginzburg-Landau equation. The tempo- 
ral Gaussian pulse shape is described by 



U(z,t T ) = U(z)exp < - 



2T 2 {z) 



(10) 



with the complex amplitude U, see Eq. (j4j. The relevant 
equations of motion are here given by Eqs. (|5bl) . (|6b[) . ([7J 
and ([8]) with B = and the nonlinearity coefficients c a = 
\/2/4, Cff, = 5^2/8, which are the same as obtained by 
the method of minimum weighted square mean error ^ 
Table 1 contains an overview of the suitable Gaussian 
test function, the relevant equations of motion and the 
coefficients for the spatiotemporal, spatial and temporal 
dynamics. 

For 5 = 0, the Gaussian is an exact solution, and 
thus the variational principle yields the exact result. For 
5 y^z 0, the Kerr nonlinearity results in an additional in- 
tensity dependent spatial and temporal chirp, see Eq. ([B]), 
and phase shift, see Eq. (jHJ). Eqs. ©-([8]) look rather com- 
plicated. However, the physics can be made rather obvi- 
ous by casting those formulas into mapping matrices for 
the complex q parameters of Gaussian bullets, see Section 
IIIII In the following, we extend the equations of motion 
to take into account arbitrary gain profiles by introduc- 
ing effective parabolic gain coefficients. In this case, the 
Gaussian is not an exact solution of Eq. ([T]) even if 5 = 0. 



B. Nonparabolic Gain Profile 

The equations of motion for the parabolic gain pro- 
file given in Eq. ^ can be modified to describe a gen- 
eral transverse and spectral gain dependence g {z, oj, x, y), 
where u is a relative frequency coordinate, centered 
around the carrier frequency ujq. The inherent symmetry 
properties of the Gaussian ansatz make it particularly 
suited for describing the propagation in media with gain 
profiles which are symmetric around x = 0, y = and 
uj = 0. Then a second order Taylor expansion of g results 
in a parabolic gain profile. However, this parabolic ap- 
proximation is only viable if the transverse and spectral 
pulse width is narrow as compared to the gain profile. In 
laser media, where a large spatial overlap with the gain is 
desired, this assumption generally fails. Also the spectral 



pulse width can significantly exceed the gain bandwidth, 
especially for few-cycle laser pulses. 

The gain term of the nonlinear Schrodinger equation, 
Eq. ([!]) , is here given by 



Q = iJT 1 {gT t {U}} , 
with the definition of the Fourier transform 

F t {U} = U = / dt r U exp (iwi r ) 



(11) 



(12) 



For general gain profiles, the Gaussian test function is 
only an approximate solution even for 5 = 0. In Ap- 
pendix A, the Kantorovitch method is used to extract the 
equations of motion for a general gain profile, which can 
be brought into the form Eqs. ©-(|S]) by defining g u , x ,y 
and go as functions of the position dependent Gaussian 
parameters: 
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(w 2 x - 2x 2 ) g 

- 2y 2 ) 9 
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1 1 
2^E 



OO /» OO 



OO J —OO 
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u 



durdxdy 



9x 2 1 9 y 2 , 9ui 2 



dujdxdy, 
(13a) 

i 

duidxdy, 
(13b) 

(13c) 
(13d) 



with the pulse energy E = n 3 / 2 A 2 Tw x w y and the spec- 
tral 1/e width defined as 



Here, 



U 



u 



n = 

is given by 
,2ttT 



y2 



Ab 2 T 2 . 



(14) 



A 



n 



exp(-x 2 /w 2 x ~y 2 /w 2 y -u; 2 /n 2 ). (15) 



Eq. (fT3"|) provides a position dependent effective parabolic 
profile, which depends on g as well as the spectral and 
transverse pulse widths. If we are interested in the purely 
spatial propagation of a Gaussian beam, Eq. (JSJ , in a 
medium with a gain profile g (z,x,y), we have to use the 
effective parabolic gain parameters 



1 



9x 
9y 
3o = 



1 



00 p co 



2x 2 )g\U\ 2 dxdy, 



— OO J —00 

OO pOO 



(w 2 y ~2y 2 )g\U\ 2 dxdy, 
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OO J —OO 
OO /"OG 



— OO J — OO 



g \Ufdxd y +^wl + ^wl, 



(16a) 
(16b) 
(16c) 



4 



TABLE I. Gaussian test function, relevant equations of motion and coefficients for spatiotemporal, spatial and temporal 
dynamics. 



Propagation Pulse Equations of Motion Coefficients 

Spatiotemporal Eq. © Eqs. ©,©,(0),© c = ^2/8, c = 7^/16 

Spatial Eq. © Eqs. ®,(Ba},(2) ) © c a = 1/4, c^, = 3/4, V = g u = 

Temporal Eq. fig) Eqs. (|5b|) . ((6b|) , dH ),© c a = y/2/4, e = 5^2/8, g = g^, = 



with the power P — ttA w x w y and 

\U\ 2 = A 2 e W (-x 2 /w 2 x -y 2 /w 2 y ) 



(17) 



For the purely temporal dynamics of a Gaussian pulse, 
Eq. (|10[) , the effective parabolic gain parameters for a 
gain profile <? (z,u;) are given by 



.9. 



1 1 



(n 2 - 2uj 2 ) g 



u 



du>, 



1 1 



U 



du + ^-n 2 , 



with the fluence F = ir 1 / 2 A 2 T and 



U 



A 2 



,2nT 



exp (-uj 2 /n 2 



(18a) 
(18b) 

(19) 



As an example, let's consider a Gaussian gain profile 



9 = .9 exp I 



A~ 2 y 2 A Z W) , 



(20) 



with the transverse 1/e gain widths A^ and A y and the 
1/e gain bandwidth A w . The equations for the effective 
parabolic gain profile, Eq. (IT5|) . evaluate to 



5 ,=5A- 2 (Aj 2 2 + l)- 3/2 (A- 2 wl + l) 



x(A^ + l)" /2 . 
, 9;c =.gA- 2 (Aj 2 fi 2 + l) 

x(A^ + l)- 1/2 , 



-1/2 



1) 



-1/2 



-3/2 



g y = gA- 2 (Az 2 n 2 + 1 

-3/2 



-1/2 



J/ 

x (A; 2 - 2 . 



A" 2 m; 2 



1) 



-1/2 



(21) 



1 



1 



.90 = g (5^ 2 + 2A 2 ) + 5 , (5™ 2 + 2A 2 ) 
+ <7„ (5w 2 + 2A 2 )]. 
With the transverse gain profile 

5 = 3exp(-A~ 2 ;r 2 -A-V), 
Eq. (|16p for the spatial beam propagation gives 



(22) 



g x = gA- 2 (A" 2 ™ 2 + l) 3/2 (A" 2 ™ 2 + l)" 
g y = gA^ 2 (A- 2 w 2 x + ly 1 ' 2 (A^w 2 , + l)^ 2 . 
go = g x (w 2 x + A 2 /2) + g y ( W 2 + A 2 /2) , 



,-i/2 



(23) 



and with the spectral gain profile 

g = gexp(~AZ 2 u J 2 ), (24) 
Eq. (|18[) for the temporal pulse propagation becomes 



g u = gA- 2 (AZ 2 n 2 + 1 



go = g u [ -W + At 



(25) 



The implementation of a general gain profile opens up 
the possibility to include a broad range of saturation ef- 
fects. For example, the equations of motion can be cou- 
pled to a differential equation for the gain profile g, de- 
scribing the evolution of g in dependence on the pulse 
parameters. 



III. DESCRIPTION BY OPTICAL MATRICES 

It is convenient to recast the equations of motion in 
a form consistent with the familiar q parameter anal- 
ysis for Gaussian optical propagation, where the opti- 
cal elements are described by matrices j 15 i 16 ' 22 This no- 
tation is compact and very practical, since it allows for 
a straightforward treatment of the successive propaga- 
tion through linear optical elements and nonlinear Kerr 
media. Additional effects in the Kerr medium, like a 
parabolic refractive index profile, can be easily incorpo- 
rated by adding the corresponding matrix elements. On 
the other hand, by rewriting the variational equations as 
differential equations for the q parameters, ABCD ma- 
trices for the Kerr effect and a nonparabolic gain pro- 
file can be extracted. We note that the matrices for 
the spatiotemporal Kerr effect obtained here are different 
from the ones derived previously based on a Taylor ex- 
pansion approach, which significantly overestimates the 
nonlinearity^ 

Using the complex q parameters q x , q y and q t , we can 
write the spatiotemporal Gaussian ansatz as 



U(z, t T , x,y) = U (z) exp 



ik y 2 



iuot 2 



_2q x (z) 2q y (z) 2q t (z) 
Comparison with Eq. ([3]) yields 



(26) 



1 



2a r 



(27) 



5 



with p = x,y, and 



<h 



1 

U} 



1 

2^2 



2b 



(28) 



Here, and q y are the reduced q parameters, with the 
vacuum wavenumber fco m the exponent of Eq. (|26l) , in- 
stead of the wavenumber in the medium. This definition 
has the advantage that a z dependent refractive index 
does not lead to additional terms in the equations of mo- 
tion for the q parameter formalism. Likewise, qt is the 
temporal analogon to the spatial reduced q parameters,— 
with o>o in the exponent, instead of the dispersion coef- 
ficient. Here, the dispersion is allowed to depend on the 
position z, without modifications of the equations. 

With the q parameters, the Gaussian beam profile for 
purely spatial propagation, Eq. ©, can be written as 



U = t/exp 



ik x 2 ik Q y 2 



2q x 2q y 



(29) 



and the temporal Gaussian pulse for describing the 
purely teporal dynamics, Eq. (fTDj) . becomes 



U = C/exp 



2q t 



(30) 



In Appendix B it is shown that the propagation equations 
for the q parameters and the amplitude can be written 
as coupled differential equations 



d z q s 

with s = x, y, t, and 



-£c' s + B' s 



2q t 2q x 2q y 



(31) 



(32) 



The coefficients B' s and C' s , which in general depend 
on the position z, can be interpreted as elements in an 
ABCD matrix of the form 



M, = 



1 B' 8 5z 
CM l 



(33) 



describing the propagation in the medium through an 
infinitely small section with length Sz. In a gain medium 
with Kerr nonlinearity, they are given by 



1 



BL = 2BU = — 



with p — 

i-woStUr 1 }) 



3 



C 



x,y, 

-1/2 



2Vu 
2igp _ 
k 

2cJ_ 
u T 2 



"o 
- 2iw 5 w , 
2c a S 



k w 2 



U 



u 



(34) 



-1/2 



w p = (fcoSjgp 1 }) T 
and the nonlinearity coefficients 



listed in Table 1. For the complex on- axis transmission 
coefficient, we obtain 



a = g + ict/,5 



U 



(35) 



The purely spatial propagation equations for a Gaussian 
beam, Eq. (|2"9"]l . are given by Eq. (|3"Tj) with s — x,y and 
Eq. (|32|) with B' t — 0. The evolution of a purely temporal 
Gaussian pulse, Eq. ([30| . is described by Eq. ([31]) with 



s = t and Eq. 032) with B' x = B' = 0. 



TABLE II. Optical matrix elements for spatiotemporal Gaus- 
sian pulse propagation through Kerr media with a spatial and 
spectral gain profile. B' t — C' t = for purely spatial beam 
propagation; likewise, B' p — C' p = for purely temporal pulse 
propagation. The Kerr parameters c a and c,/, are listed in 
Table 1. A general nonparabolic gain profile can be taken 
into account by using the g^, g v and go given in Eqs. (|13|) . 
(|16[) and (JTHJ) for spatiotemporal, purely spatial and purely 
temporal propagation, respectively. 



Effect 


B p 


Cp 


B[ 


c\ 


a 


Free space 


no 1 












Dispersion 






2Vlo 








Parabolic gain 





2ig p k 


2ig,jLo 





go 


Kerr effect 










C a - L 


j^\U\ 2 ic,6\U\* 



The coefficients are composed of various contributions, 
B' s = J2i B' s i , C' s = J2i Csf where the index i represents 
the different physical mechanisms affecting the propaga- 
tion. Each of these effects itself can be described by a ma- 
trix of the form Eq. (|33|) . We identify the matrix elements 
for free space propagation (B p .i = Uq 1 5z and C v .\ = 0), 
soft aperturing (B P 2 — and C p ,2 = 2ig p /c ~ 1 <5z), dis- 
persion (B tt i = 2T>loq5z and C^i = 0), and spectral 
parabolic filtering (B tj2 = 2iu} g UJ 5z and C t ,i = 0)^^^ 
The Kerr effect is incorporated as an intensity dependent 
lens2£ in the spatial domain, describing the self-focusing 
action, and an intensity dependent "chirper" or temporal 
lens2£ in the time domain for the self-phase modulation. 
An overview of the matrix elements for the different ef- 
fects is given in Table 2. Note that a general gain profile 
deviating from the parabolic approximation can be in- 
corporated by using above gain matrix elements together 
with the equations for g p , g^ and go given in Section lll Bl 

Additional effects can easily be incorporated by adding 
further matrix elemens. For instance, a parabolic refrac- 
tive index profile n (r) = no — {fi2^ x x 2 + U2 :V y 2 ) /2, as 
generated by thermal lensing in a laser rod, can be taken 
into account by the elements B' p = 0, C' p — —ri2, p - Also 
transversely varying saturable gain can be treated by ad- 
ditional ABCD matrices^ 



IV. EXAMPLE 

An important application of the Gaussian approxi- 
mation described above is the simulation of the opti- 
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cal propagation in laser resonators. Examples are the 
temporal pulse evolution in mode-locked fiber lasers, the 
laser beam propagation in high-power laser rods, or the 
spatiotemporal pulse dynamics in Kerr-lens mode-locked 
lasers. Here, the steady-state solutions cannot be ob- 
tained directly from the optical matrices, since the ma- 
trix elements for the Kerr nonlinearity and the general 
gain profile depend on the pulse parameters. Rather, 
the steady state must be found by iteratively solving the 
equations, i.e., propagation over many roundtrips till the 
steady state is reached. 




L, I, L l 2 L 2 

FIG. 1. Simplified model of a Kerr-lens mode-locked laser 
resonator. The end mirrors are represented by dashed lines. 



As an example, we choose here the spatiotemporal 
pulse propagation in a Kerr-lens mode-locked laser res- 
onator. The setup, shown in Fig. 1, consists of a non- 
linear Kerr medium and linear resonator arms to its 
left and right, which contain an element with nega- 
tive dispersion and a focusing element.— For the res- 
onator arms, we assume a focal length /z, = 5 cm and 
a group delay dispersion Si — S2 = — 75 fs 2 . For the 
Kerr medium, the material parameters of Tksapphire are 
used, with V = 60fs 2 /mm, 5 = k n 2 i = 0.25 pm/ MW, 
n = 1.76, f = 375 THz and B = l/(2n fc ) = 36.2 nm. 
The lengths in the resonator are given by L\ = 80 cm. 
L2 = 110 cm, li = 5.05 cm, L = 0.25 cm and I2 = 5.2 cm. 

In Ref. 4, the Gaussian solutions for the setup in 
Fig. 1 are obtained taking into account only the energy- 
conserving effects, i.e., neglecting gain and loss. Here, 
a Gaussian gain profile in the Kerr medium is added, 
characterized by the parameters g, A w and = A y , 
see Eq. (|2TJ|) . The output coupling is taken into account 
by normalizing the intracavity pulse energy at the right 
end mirror to E = 20 nJ. In the following, we examine 
the dependence of the solution on the gain parameters. 
Fig. 2 shows the Gaussian pulse duration T at the right 
end mirror and the stability factor T as a function of 
the roundtrip gain G = AE/E for different values of A x 
and A u . T is defined as the ratio between the gain per 
roundtrip in pulsed and cw operation, T — G/G cw , where 
G cw is the roundtrip gain for the Gaussian steady-state 
beam solution of the linear resonator. The stability fac- 
tor serves as an indicator for the suppression of the cw 
lasing in pulsed operation, a value of T > 1 indicating 
stable pulsed operation^ 

Fig. 2(a) shows the pulse duration T at the right end 
mirror. The pulse duration increases for strong spectral 
gain filtering, i.e., for a strong roundtrip gain G in com- 
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FIG. 2. Steady-state pulse solution of the setup shown in 
Fig. 1 for a fixed intracavity pulse energy of 20 nJ at the 
right end mirror. Shown is a) the pulse duration T at right 
end mirror and b) the stability factor T as a function of the 
gain per roundtrip G. Results for Aw/ (2ty) = 40, 45, 56 THz 
and A w — > 00 are represented by dotted, dash-dotted, dashed 
and solid lines, respectively. Transverse gain widths of A x — 
10 /im, A x — 20 /jm and A x — > 00 are marked by light gray, 
dark gray and black colors. 

bination with small gain bandwidth A w . Also a small 
transverse gain width A x leads for a fixed value of G to in- 
creased spectral filtering, mainly because of the increased 
peak gain g. Thus, minimum pulse durations can be ob- 
tained in a setup with a small gain and weak outcoupling 
at the end mirrors. Fig. 2(b) shows the stability param- 
eter r. Due to the high peak power, the pulsed solution 
experiences self-focusing and thus transverse contraction 
in the Kerr medium, leading to an increased transverse 
overlap with the gain. On the other hand, the cw solution 
is not affected by spectral gain filtering. Thus, T reaches 
maximum values for a small transverse gain width A x 
and a broad gain bandwidth A w . 

The results for vanishing gain, G — > 0, coincide with 
the ones obtained in Ref. 4, taking into account only 
the energy-conserving effects. For a small gain of a few 
percent, the energy-conserving dynamics is still a good 
approximation, as can be seen from Fig. 2. However, the 
equations of motion with gain have the distinct advantage 
that the system is attracted by its stable solutions, while 
for the energy-conserving equations, additional boundary 
conditions have to be introduced to find the steady-state 
solutions^ 

As pointed out in Ref. 4, higher-order dispersive ef- 
fects, which are not included here, can considerably affect 
the pulse shape. In addition, the Gaussian approxima- 
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tion fails for excessive self-focusing and self-phase mod- 
ulation, as well as strong nonparabolic gain aperturing 
and filtering. Under extreme conditions, the nonlinear 
Schrodinger equation itself, as given in Eq. ([1]) , loses its 
validity. 



V. CONCLUSION 

In conclusion, we have studied the spatial, temporal 
and spatiotemporal optical propagation in Kerr media 
with transverse and spectral gain filtering by applying 
the variational principle. Based on the Kantorovitch 
method, we derived the Gaussian equations of motion 
for parabolic and general gain profiles. By reformulat- 
ing the variational equations as differential equations for 
the q parameters, we could extract ABCD matrices for 
the Kerr effect and a general transverse and spectral 
gain profile. As an example, we studied the steady-state 
spatiotemporal Gaussian pulse dynamics in a Kerr-lens 
mode-locked laser resonator. 

The equations of motion can be solved efficiently with a 
standard differential equation solver and allow for a quick 
simulation of the Gaussian optical propagation through 
gain media with Kerr nonlinearity. By iterative solution 
of these equations, the steady-state pulse or beam shape 
in a laser resonator can be obtained. Further effects, 
like a parabolic refractive index profile as generated by 
thermal lensing in a laser rod, can easily be considered 
by additional ABCD matrices. Gain saturation can be 
taken into account by complementing the equations of 
motion with suitable gain saturation equations. 
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APPENDIX A: DERIVATION OF THE VARIATIONAL 
EQUATIONS 

In this appendix, we derive from Eq. ([1]) the equations 
of motion for the spatiotemporal Gaussian pulse parame- 
ters, using the Kantorovitch method. 12 The conservative 
Lagrangian is given by 9 



C 



U 



dU 



u 



dU* 



(Al) 



while the non-conservative process is described by the 
expression Q, given in Eq. (j5J| and Eq. respectively. 
For the envelope U, we insert the test function Eq. ([3]). 



The Euler-Lagrange equations for the real parameter 
functions f — A,<p, T, w Xl w y , 6, a x , a y are then given byi 2 . 





dU 
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dU 


2 


dU 




-B 


-B 


dt r 


dx 


dy 



d(C) d d(C) 

T = H f> 



df dz df 
with the reduced Lagrangian 



CO 



oo poo poo 



— oo J — oo J — oo 



£dt T dxdy 



and the non-conservative term 

pOO /*00 /■CO QJJ 



Rf = 23? 



oo J — oo J — oo 



(A2) 



(A3) 



(A4) 



Using the definition of the Fourier transform in Eq. (JT2J) 
and Parseval's theorem, we can with J- t {Q} — igU ex- 
press Eq. (|A4[> as 



1 If 00 f°° f°° ~ dU* 
R S = Z^VI I I gU-Qjrdudxdy \ . (A5) 



— oo o — oo •/ — OO 



Since g is real, we have Ra = Rw x = Rw y — 0. Fur- 
thermore, for a parabolic gain profile with the gain term 
Eq. ([2j , we obtain 

R<t> = E [2g - g x w 2 x - g y w 2 y - gM 2 } , (A6) 



R T =4g u EbT-\ (A7) 



R b = -ET 2 [2g - g x w 2 x - g y w 2 y + g u (T" 2 - 12T 2 6 2 )] , 

(A8) 

Ra x = \ew 2 x [2.9o - 3g x w 2 x - g y w 2 y ~ gM 2 } , (A9) 

and a corresponding expression for R ay , where the pulse 
energy is given by 



E = ir 3/2 A 2 Tw x w y . (A10) 

From Eqs. (|A6|) - (IA9|) . we can extract equations for g^, 
g x , g y and g , 

gu = RtT / (AEb) , (All) 
g x = (R c ,w 2 x /2-R ax )/(Ewt), (A12) 



g y = (R^w 2 y /2-R ay )/(Ew 4 y ), 



(A13) 



.90 = (Rt/E + g x w 2 x + g y w 2 y + gM 2 ) /2. (A14) 

This enables us to express any gain profile formally 
through effective parabolic gain parameters, which de- 
pend on both the gain function and the pulse parameters. 
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Inserting Eq. (|A5|) into Eqs. (|A11[) - (|A14I) . we arrive at 
Eq.(S). 

Setting / = j4 in Eq. jA2]) yields 



and for / = 



2<\>' - (6T 2 + oX + a^) + 2? + 4fo 2 T 2 ^ 

2„„2 N \ _ g / 1 j_ /l^,^ , ;A 



4^ 



*7I = > 



(A15) 



E' = E 



we get 



25o ~ 9xw 2 x -g y w I- g 



y2 



46 2 T 2 



(A16) 

The Euler-Lagrange equations for the pulse duration and 
the beam widths are given by 



APPENDIX B: DERIVATION OF THE EQUATIONS FOR 
THE Q PARAMETERS 

The spatiotemporal dynamics is described by coupled 
equations of motion for q Xl q yi q t and U. Differentiation 
of Eq. (|27J) with respect to z yields with Eqs. (|5a| and 
(|6a[) the differential equation for q p , 



2i#p 
k 



2c a S 
k wl 



U 



1 

n ' 



(Bl) 



with p — x,y and w p — (fco^ {^p -1 }) • The differen 
tial equation for q t is obtained by differentiating Eq 
with respect to z and inserting Eqs. (|5b|) and (|6bl : 



b'T 2 - 3a' x w 2 x - a' y wl +V(^+ Ab 2 T 2 ) 



-B ( + Walwl ) - fi f J, + 4a 2 w 2 



(A17) 



b'T 2 - a' x w 2 x - Za' y w 2 y + V + 46 2 T 2 ) 



I2a 2 y w 2 y 



6 A2 _ n 

(A18) 



- 36'T 2 - a' x w 2 x - a' y w 2 y + V (- ^ + 126 2 T 2 ) 



-B f J* + 4a 2 x w 2 x ) - B ( A 



4a 2 W 2 



and for the chirp parameters, we obtain 

1 w' 

E ' + E-*- - ABEa x 

2 W r 



(A19) 



E 



go 



(A20) 



1 , w' 

-E 1 + E-± - ABEa v 

2 w„ 



= E 



gx 2 3g y 2 
9o - yw x - — w y 



(A21) 



-E' + E— + AVEb 
2 T 



E 



9x 2 9y 2 , n 

9o ~ -ttW x - — w a + 2g u 



(455 - 3T262 



(A22) 



Inserting Eq. (|A"l6|) into Eqs. (|A"20|) - (TA"2"2"[) results in 
Eqs. pp and ([5b}. Multiplying Eqs. ([AlTl) . ([A18|) and 
(|A19|) by 2 and subtracting Eq. (|A15|) from each equa- 
tion yields Eqs. ([6a|) and ([6b)) . Eq. (|7|) can be obtained 
from Eq. (|A16[) by inserting Eqs. (|A10I) . ((5aj) and ([5b]l. 
Multiplying Eq. (|A15I) by a factor of 2^ and subtracting 
Eqs. (|AT7|) . (|AT8t> and ([ATI) yields Eq. ©. 



3< 



2c a <5 
^oT 2 



[7 



q 2 + 2Vlo + 2i 5w wo, (B2) 



with T = (— cjoSj'Zt *}) 



-1/2 



Differentiating Eq. dU 



with respect to 2 and inserting Eqs. (JJJ and (jHJ yields 
with Eqs. (|2T|) and (|28|l the equation of motion for the 
complex amplitude, 



U A 



Bk Bk Vuj ,ui g u 



Qt 



-go+ic^S 



U 



(B3) 

The purely spatial beam propagation, Eq. (|29|) . is de- 
scribed by Eqs. ([Bl]) and jB3j with = ^=0, and 
the temporal pulse propagation, Eq. (15U1) , is described by 
Eqs. ([El) and JB3]) with B = 0. The nonlinearity coeffi- 
cients are given in Table 1. 

In the q parameter formalism, discrete optical elements 
are represented by 2 x 2 matrices 



Mn = 



B s 



(B4) 



For the spatiotemporal pulse propagation, each optical 
element is characterized by three matrices, i.e., s — 
x,y,t. The transformation law for the propagation 
through an optical element extending from position z\ 

to Z2, 



9s (22) 



A s q s (zi) + B s 
C s q s (zi) + D s 



(B5) 



is valid for both the spatial and temporal q 
parameters . 15 i 16 The amplitude at position z-i is given 
by 

U (z 2 ) = tU (zi) [A x + B x /q x ( Zl )]- 1/2 

X [A y + B y /q y ( Zl )}- 1/2 [A t + B t /q t ( Zl )r 1/2 

(B6) 

for a spatiotemporal Gaussian pulse, 

U (z 2 ) = tU (zi) [A x + B x /q x ( Zl )]- 1/2 

x[A y + B y /q y ( Zl )r 1/2 (B7) 
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for a spatial beam, and 

U (z 2 ) = tU (zi) [A t + B t /q t ( Zl )}- 1/2 



(B8) 



for a purely temporal pulse, with the on-axis trans- 
mission t = exp(J'adz), where a is the complex on- 
axis transmission coefficient. The propagation equations 
Eqs. (|B1|) - (|B3[) can be obtained by dividing the Kerr 
medium into small sections of length Sz, and represent- 
ing each section by ABCD matrices of the form Eq. (f33|). 
From Eqs. (|B5|) and (|B6[) . we obtain 



q s (z + 5z) 



and 



q s (z) + B' s Sz 
q s C' s Sz + 1 



-q 2 s (z) C' s Sz+q s (z)+B' s Sz 
(B9) 



U (z + Sz) = U (z) exp (a5z) (1 + B'Jz/q x ) 1/2 

x (1 + B' y Sz/q y y 1/2 (1 + B' t 8z/q t )- 1/2 

f Tf x (, B'Jz B' y Sz B' t 8z 

~U(z)[l- -f ^— + aSz 

\ 2q x 2q y 2q t 

(BIO) 

In the limit Sz -t 0, this results in Eqs. d3U) and (1521) . 
Comparison with Eqs. (iBTj). (|B2l) . and (|B3|> yields the 
elements Eq. (|34|) and the a given in Eq. (|35l) . The equa- 
tions for the purely spatial or temporal dynamics can be 
derived in an analogous manner. 
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